arXiv:1506.06478vl [gr-qc] 22Jun2015 


Slowly rotating regular black holes with a charged thin shell 
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We obtain rotating solutions of regular black holes which are constructed of de Sitter spacetime 
with the axisymmetric stationary perturbation within the timelike charged thin shell and the Kerr- 
Newman geometry with sufficiently small rotation outside the shell. To treat the slowly rotating 
thin shell, we employ the method developed by de la Cruz and Israel. The thin shell is assumed to 
be composed of a dust in the zero-rotation limit and located inside the inner horizon of the black 
hole solution. We expand the perturbation in powers of the rotation parameter of the Kerr-Newman 
metric up to the second order. It is found that with the present treatment, the stress tensor of the 
thin shell in general has anisotropic pressure, i.e., the thin shell cannot be composed of a dust if the 
rotational effects are taken into account. However, the thin shell can be composed of a perfect fluid 
with isotropic pressure if the degrees of freedom appearing in the physically acceptable matching of 
the two distinct spacetimes are suitably used. We numerically investigate the rotational effects on 
the spherically symmetric charged regular black hole obtained by Uchikata, Yoshida and Futamase 
in detail. 

PACS numbers: 04.70.-s 


I. INTRODUCTION 

The interior of the black hole is still an open problem in 
gravitational physics. Although the singularity theorems 
[Ml predict the inevitable formation of the spacetime 
singularity due to the gravitational collapse to the black 
hole, the theorems are valid only if certain conditions 
hold within the classical framework. For example, the 
strong energy condition is assumed to prove the singu¬ 
larity theorems. This condition is given by R^RP > 0, 
where R and U * l are the Ricci tensor and a nonspace¬ 
like vector, respectively. In de Sitter spacetime, which is 
a vacuum solution of Einstein equations with a positive 
cosmological constant, A > 0, the strong energy condi¬ 
tion is violated since R pv Pl v = A < 0 for Pl^ < 0. 
However, the so-called weak energy condition holds in 
this spacetime. Therefore, we could construct black holes 
without spacetime singularities, i.e., regular black holes 
or nonsingular black holes, using spacetime with a posi¬ 
tive cosmological constant. 

In this study, we focus on regular black hole models 
containing a regular center inside the black hole event 
horizon [4[. Regular black holes of this type basically 
have causal structures similar to that of the Reissner- 
Nordstrom black hole but with regular centers instead 
of spacetime singularities a®. Near the center, they 
behave like the de Sitter spacetime. Similarly to the 
de Sitter spacetime, in these regular black holes, the 
weak energy condition is satisfied everywhere, but the 
strong energy condition is violated near the center. Thus, 
these regular black hole models are outside the scope of 
the singularity theorems. In other words, the spacetime 


* nami.uchikata@ist.utl.pt 

I yoshida@astr.tohoku.ac.jp 


singularity-free property of these regular black holes is 
not inconsistent with the singularity theorems. 

The models of the regular black hole we consider can be 
divided into two classes from the aspect of smoothness of 
their spacetime. The models belonging to one class are 
characterized by a smooth metric, which was first sug¬ 
gested by Bardeen [§]. In the present paper, we call the 
regular black hole of this class the Bardeen-type regu¬ 
lar black hole. The solution is typically described by the 
gravitational field coupled to the nonlinear magnetic field 
[MlH- A magnetic monopole behaves like the positive 
cosmological constant so that the spacetime can become 
regular everywhere if there is a magnetic monopole inside 
the event horizon near the center of the solution. A key 
ingredient for the solutions of this class is the existence 
of fields acting as an effective positive cosmological con¬ 
stant inside the event horizon. This class of regular black 
hole solutions was generalized to the unc harg ed cases 
flTO , charged case [15] and rotating case |l6l ] Forma¬ 
tion and evaporation of Bardeen-type regular black holes 
were considered by Hayward flil ]. Quasinormal modes 
of Bardeen-type regular black holes were investigated in 
Ref. Q3. 

The regular black hole model belonging to the other 
class is characterized by its spacetime composed of two 
distinct spacetimes: one is de Sitter-like spacetime and 
the other is black hole spacetime, matched with a thin 
shell inside the event horizon. Models of this class are 
motivated by the conjecture suggesting that the space- 
time curvature_ has an upper limit of the order of the 
Planck scale [H|-[20|; e.g., the curvature invariant is re¬ 
stricted by Rfj.i/a/3R^ l ' a ^ < l ~ 4 with R^a/3 and l p being 
the Riemann tensor and the Planck length, respectively. 
According to this conjecture, the collapsing matter turns 
into a de Sitter phase when the curvature reaches the up¬ 
per limit. For the regular black hole model of this class, 
as mentioned before, de Sitter spacetime, whose curva- 
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ture invariant is given by R v , V a 0 R^ va ^ = 24/L 4 with L 
being the de Sitter horizon radius, is frequently consid¬ 
ered as a regular core instead of a spacetime singularity. 
This conjecture suggests that the de Sitter horizon ra¬ 
dius L for this model is of the order of the Planck length 

mm. 

In this study, we focus on the regular black hole of the 
second class. For the regular black hole solutions of this 
second class, the uncharged I2ltj23| and charged spheri¬ 
cally symmetric solutions [24|,[25| have been investigated. 
The radial stability of these solutions is studied and it 
is shown that some regular black hole solutions are sta¬ 
ble against small radial disturbances [H, |25[ . Although 
many studies on the regular black hole of the second class 
have so far focused on spherically symmetric nonrotating 
cases, it is reasonable to generalize the regular black hole 
to rotating cases because astrophysical objects usually 
have angular momentum due to its conservation prop¬ 
erties. In this study, therefore, we will extend previous 
analyses of spherically symmetric regular black holes to 
rotating cases, i.e., cases of the rotating regular black 
hole of the second class. 

As a first step, we focus on the slowly rotating case. 
Rotational effects are then treated as small perturbations 
around spherical solutions. More specifically speaking, 
we consider the case of |a| -C M, where a and M are the 
rotation parameter and mass of the black hole solution. 
The rotational effects up to the second order of \a/M\ 
are taken into account in this study because we are con¬ 
cerned with rotational effects on the matter properties of 
the thin shell like deformation of the shell. We assume 
the thin shell to be a timelike hypersurface as done in 
Ref. (2f|. In order for the timelike stationary thin shell 
to exist inside the event horizon, we need the inner hori¬ 
zon besides the event horizon because such a timelike sta¬ 
tionary thin shell has to be inside the inner horizon. This 
implies that a charged black hole has to be considered as 
the solution outside the shell because a slowly rotating 
Kerr black hole within the accuracy we assume has the 
event horizon only. In this study, therefore, we assume 
the charged regular black hole obtained in Ref. [25[ to be 
the unperturbed spherically symmetric solution. We ba¬ 
sically use a method of solution similar to that given by 
de la Cruz and Israel [26j for matching between the Kerr 
spacetime and the regular vacuum solution of Einstein 
equations under the assumption of slow rotation. How¬ 
ever, not only the spacetime but also the electromagnetic 
fields have to be considered in the present case because 
the thin shell is charged. Outside and inside the shell, 
we assume a slowly rotating Kerr-Newman solution and 
a stationarily and axisymmetrically perturbed de Sitter 
solution, respectively. The perturbed quantities associ¬ 
ated with the thin shell are derived from the junction 
conditions for the spacetime and electromagnetic field 
(see, e.g., Refs. 1271 - 1291 ). By using analytically obtained 
results, we show some numerical examples for the slowly 
rotating charged regular black holes. 

The plan of this paper is the following. : In Sec. II, 


we give the formalism for obtaining the slowly rotating 
charged regular black hole and the analytical results de¬ 
rived in the present study. In Sec. Ill, some numerical 
results for the slowly rotating charged regular black hole 
are shown. A conclusion is given in Sec. IV. Throughout 
this paper, we use the units of c = G = 1, where c and 
G are the speed of light and the gravitational constant, 
respectively. 


A. 


II. METHOD OF SOLUTION 

Basic equations for matching of two distinct 
spacetimes 


As mentioned before, we construct single spacetime in 
which two different spacetimes are matched with a thin 
shell. For the reader’s convenience, here, we concisely 
describe the formalism developed in Ref. [13, H3 > which 
gives a coordinate-independent prescription for match¬ 
ing of two different spacetimes. Let V ± be the four- 
dimensional spacetimes that have metrics g^ and sys¬ 
tems of coordinates x ±a . Let E be a three-dimensional 
timelike hypersurface described by intrinsic coordinates 
y a and located at the boundaries of V + and V~. Thus, 
E is given by x ±a = x ±a (y a ). Here and henceforth, we 
use greek and roman lower case letters to describe in¬ 
dices of the four-dimensional spacetime and of the three- 
dimensional hypersurface, respectively. Let n a be the 
unit normal vector to the timelike hypersurface E. Thus, 
n a has to satisfy 

n a n a = 1, e“n a = 0, (2.1) 

where e“ is the basis vector on E, defined by 

_ ux 
a dy a ' 

Here and henceforth, ‘±’ is sometimes omitted for the 
sake of brevity. The induced metric h ab and the extrinsic 
curvature K a b associated with E are, respectively, defined 
by 

hab = ffa/?e“ef , K ob = -n a;j ae“ef. (2.3) 

Here and henceforth, we denote the covariant differenti¬ 
ation associated with g a p and h ab by the semicolon (;) 
and the stroke (|), respectively. So that V + and V~ are 
smoothly joined across E, the first junction condition, 
given by 


( 2 . 2 ) 


[hab] = 0 , 


(2.4) 


where [g] = g| + — q\~, and q]^ is the value of the physical 
quantity q evaluated on E by taking the limit from V±, 
has to be fulfilled. Jump in the extrinsic curvature on E 
relates to the energy-momentum tensor on E, S a b , which 
is explicitly given by 

Sab = - 5 — {[Kab] ~ h a b[K ]) . 

57r 


(2.5) 
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When S a b ^ 0, we regard E as a thin shell because there 
is the matter field on E. Equation E3 is called the 
second junction condition in this paper. 

When there are electromagnetic fields in , the junc¬ 
tion conditions for electromagnetic fields are also required 
so that the electromagnetic fields are smoothly joined 
across E, as argued in Ref. [29|]. The junction conditions 
for the electromagnetic field are summarized as follows: 


Here, F a p stands for the Faraday tensor, given in terms 
of the vector potential A M by F a p = d a Ap — dpA a , and 
j a is the current vector on E. The first and second con¬ 
ditions in Eq. (ESI), respectively imply that the tangent 
components of the Faraday tensor to E have to be contin¬ 
uous through E and that the current on E is given by the 
jump of F an on E. These junction conditions are exactly 
the same as those considered in classical electrodynamics 
(see, e.g., Ref. [30[ )• 


B. Two spacetimes matched with a thin shell 

In order to obtain rotating regular black hole solu¬ 
tions, we match two different spacetimes by using a time¬ 
like thin shell. Inside and outside the thin shell, sta¬ 
tionary axisymmetric regular spacetime with a positive 
cosmological constant V~ and the Kerr-Newman space- 
time V + , respectively, are assumed in this study. The 
matching surface E lies inside the inner horizon in V + 
and inside the de Sitter horizon in V~ as considered in 
Ref. [25| . The effects of rotation are assumed to be suf¬ 
ficiently small that they can be treated as perturbations 
around a nonrotating spherically symmetric solution con¬ 
sidered in Ref. [25|. Here, the slow rotation means that 
the rotation parameter of the black hole is sufficiently 
small. Thus, it is useful to introduce a smallness param¬ 
eter representing the strength of the rotational effect as 
e 2 = a 2 /M 2 <C 1, where a and M are the rotation param¬ 
eter and mass of the black hole. It is physically reason¬ 
able to expect that the rotation parameter of the black 
hole solution is of the same order as the angular veloci¬ 
ties of the thin shell and frame dragging of the spacetime 
inside the thin shell. This expectation is justified when 
fully consistent solutions are obtained by the matching of 
two distinct spacetimes with a slowly rotating thin shell. 
In this study, we assume that rotating regular black hole 
solutions have equatorial symmetry and time-azimuth re¬ 
flection symmetry because of the symmetry of the Kerr- 
Newman solution. Therefore, the deformation of the thin 
shell appears as the second-order effects of e, which cor¬ 
respond to centrifugal effects due to the thin shell’s spin. 
We then obtain all the physical quantities of the regular- 
black hole within an accuracy up to the e 2 -order. 


Since rotating charged thin shells are considered, as 
mentioned before, we assume the exterior spacetime of 
the shell, V + , to be a rotating charged black hole solution 
i.e., the Kerr-Newman solution in the Boyer-Lindquist 
coordinate and expand it up to the e 2 order as follows: 


ds 2 = \ -f RN (r) + , 




cr 

II 

o 

and 

e 

II 

(2.6) 

+ 

where 




+ 

F ab = E a/3 e“ef, 

and 

F _ F aP a n 0 

z an z ap c a 

(2.7) 

+ 


/. 


1 


r 


[ Jrn{ 

+ e 2 M 2 cos 2 9)d9 
e 2 M 2 ( 



2 r 2 [ 


2 r(M + r) — Q 1 


+ (Q 2 — 2 Mr) cos 29 


2 a 


Ql , 2 2M % in 2 6)d4>dt + 0(e 3 ), 


sin Qdcjr 


( 2 . 8 ) 


where /iiAr(r) = 1 — 2 M/r + Q 2 /f 2 and Q is the charge 
of the black hole. In this study, we assume that the 
gravitational mass M for the Kerr-Newman solution can 
change due to rotational effects of the thin shell. This 
change in M is a second-order effect with respect to the 
spin parameter a. Then, we explicitly assume that M —> 
M + e 2 SM if the thin shell rotates. On the other hand, 
we require charge conservation for the thin shell i.e., Q 
is independent of the spin parameter a. This assumption 
will be explored in Sec. II.E. The vector potential, H+, 
is similarly given by 



We assume the interior spacetime of the shell, V ~, 
to be a stationarily and axisymmetrically perturbed de 
Sitter solution given in the so-called static coordinates. 
The perturbation is too expanded around the de Sitter 
solution written in the polar coordinates (f, r, 9, 4>) up to 
the e 2 order as follows: 


ds 2 = - fdsir) (l + 2e 2 h(r , 9)) di 
2e 2 m(r, 9) 


+ r 2 (l + 2e 2 k(r,9)) 

x | d9 2 + sin 2 9 ( d(f> — 0(r)df) 2 j + 0(e 3 ), 


( 2 . 10 ) 


where fds(r) = l — r 2 /L 2 with L being the de Sitter hori¬ 
zon radius and 0(r) = O(e). The functions appearing 
in Eq. (12.101) . 0(r), h(r,9), m(r,9) and k(r,9) are reg¬ 
ular solutions of perturbed Einstein equations coupled 
to a positive cosmological constant and electromagnetic 
fields generated by the charge of the shell. In order to 
achieve separation of variables and complete matching up 
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to the e 2 order between the two spacetimes described by 
Eqs. m and (12.101) . the functions h(r,9), m(r,0 ) and 
k(r,9) are expanded as follows: 

h(r,9) = h 0 (r) + h 2 (r)P 2 (cos 9), 

m(r,9) = mo(r) + m 2 (r)P 2 (cos 9), (2.11) 

k(r,9) = k 2 (r)P 2 (cos 9). 


where Pi denotes the Legendre polynomial of degree l 
(for this treatment of the metric perturbation, see, e.g., 
Refs. [Hl-dil). The vector potential A~ is assumed to be 


A~ = 


-^ + e 2 { B 0 (r ) + B 2 (r)P 2 (cos9)} 
+ aA 3 (r ) sin 2 9d(f) + 0(e 3 ), 


dt 


( 2 . 12 ) 


where R denotes the radius of the thin shell in the limit of 
e —> 0. Analytical expressions for the functions appear¬ 
ing in the solution for V~ are derived from perturbed 
Einstein-Maxwell equations with a positive cosmological 
constant in Appendix A. 


C. Equations of the thin shell S 


Here x M and aA are the four-dimensional coordinates out¬ 
side and inside the shell, respectively. The functions £(0) 
and C( 0 ) can be chosen not to violate the first junction 
condition dm as discussed later. We further assume 
that 

A = A o + e 2 A 2 (0) + O(e 3 ). (2.17) 

The tangent vectors to E evaluated outside and inside 
the shell, and e£, are then given by 

e£ = (l, 0 , 0 , 0 ) + 0 (e 3 ), 
g£ = (0,e 2 C,e,l,0 )+0(e 3 ), (2.18) 

= (0, 0 , 0 , 1 ) + 0(e 3 ), 

ej. = (A, 0,0,0) + 0(e 3 ), 

e e = ( e 2 A 2 ,e, e 2 £,o, 1 + e2 l%i 0) + 0(e 3 ), (2-19) 

e $ = ( 0 , 0 , 0 , 1 ) + 0 (e 3 ). 

The functions £ and £ are expanded as 
£(©) = & + 6 -P 2 (cos 0 ), 

C(0) = Co + C 2 ^ 2 (cos 0 ). ’ 


In order to match the two spacetimes, first of all, we 
need to specify the shape of the matching layer E. The 
matching layer or thin shell E is assumed to be given by 

r = R+e 2 £(0) + O(e 3 ), r = R+e 2 £(9) + 0(e 3 ) , (2.13) 

in V + and V~, respectively. Then, we have 

n +=N+(0,l,-e 2 ( § ,0) + 0(e 3 ), 

n~ =N-(0,l,-e 2 io,0)+0(e 3 ), 

where N± is determined by the condition of = 1 

and the explicit expressions are given by 


N+ 


fRN 


1 



1 - e" 


r 2 (M 2 — 2 8Mr) — M 2 f R N cos 2 9 
2r 2 f RN 


+ 0(e 3 ), 



m(r , 9) \ 
rfds ) 


+ 0(e 3 ). 


(2.15) 


We choose y a = (T, 0, $) as the intrinsic coordinates of 
the shell and assume that the shell, E, is given by 

x +l1 = , 

= ( 't,r,6,4bj , 

= (T, R + e 2 £(0), 0, d>) + 0(e 3 ), 
x~' 1 = x M , 

= (t, r, 9 , cj)) , 

= (AT, R + e 2 £(0), 0 + e 2 Z e (0), $) + 0(e 3 ). 

(2.16) 


D. First junction conditions 


With the assumptions made before, the nonzero com¬ 
ponents of the induced metrics on are given by 


h TT = -fRN(R) + e 2 |2£ 

+ 0(e 3 ), 

( 2 . 21 ) 

h+ Q = R 2 + e 2 (2 R( + M 2 cos 2 0) + 0(e 3 ) , (2.22) 

hta, = R 2 sin 2 0 


+M 


2 Q 2 ~2MR 2^,2 SM\ 

-r^ cos 0 + _ ? r ] 


M 2 


+ e 2 2R(+ — {2R 2 + (2MR~Q 2 )(1 
x sin 2 0 + 0(e 3 ), 

hr* = h+ T = a ® sin 2 0 + 0(e 3 ), 


cos 20)} 

(2.23) 

(2.24) 


i tt = ~ A o fds(R) ~ 2e 2 Ao } hf d s(R) - 


- 2 e 2 f dS A Q A 2 + A 2 R 2 0 2 sin 2 0 + 0(e 3 ), (2.25) 

Hqq = R 2 + e 2 (2 (R + 2 kR 2 + 2R 2 l%) + 0(e 3 ), (2.26) 
h TO = h OT = -£ 2 Aofds(R)A 2 ,e T + <J(e 3 ) , (2.27) 

/i = R 2 sin 2 0 


+ e 2 (2^R + 2kR 2 + 2R 2 ^^-l 0 \ sin 2 0 + 0(e 3 ), 
\ sin 0 ) 

(2.28) 

h~^ = -0A o R 2 sin 2 0 + 0(e 3 ). (2.29) 
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The e°-order condition is obtained from the first junc¬ 
tion condition, [ Htt] = 0 , which is given by 

A 0 2 = . (2.30) 

JdS 

The first junction condition [hr®] = 0 gives us the e-order 
condition, given by 

Q 2 - 2 MR 

“~~ a AqR* ’ (2.31) 

= auii , 

where we have used the relation ui = const., which is the 
regular solution for our problem as shown in Appendix 
A. Here, it should be emphasized that the two space- 
times can be matched with a slowly rotating thin 
shell with accuracy up to the order e because the first 
junction conditions are fulfilled up to the order e, under 
our assumptions. There is no degree of freedom we can 
specify up to this order. All the physical quantities are 
fully determined by the unperturbed regular black hole 
solution and the spin parameter a. 

The e 2 -order conditions required by the first junction 
conditions are summarized as follows: 


For the e-order quantities, then, we have the seven unde¬ 
termined constants, SM, £o, £ 2 , Co, C 2 , C 4 , and D- 2 . Note 
that C 4 and D 2 appear in the l = 0 regular solution 
(ho(R)) and the l = 2 regular solutions (h 2 (R), »n 2 (i2) 
and k 2 (R)), respectively (see Appendix A). Among these 
constants, on the other hand, the first junction condi¬ 
tions give us the four constraint equations (12.341) (12.371) . 
Equations (12.341) and (12.361) are constraints for SM, £ 0 , 
Co, and C 4 . Equations (12.351) and (12.371) are constraints 
for £ 2 , ( 2 , and D 2 . Therefore, we can freely choose two 
constants among SM, Co, Co, and C 4 , and one constant 
among £ 2 , C 2 , and D 2 not to violate the four constraint 
equations (12.341) (12.371) . In principle, these three degrees 
of freedom in the set of the seven constants, SM, Co, 
£ 2 , Co, C 2 , C 4 , and D 2 , which completely determine the 
structure of the slowly rotating regular black hole con¬ 
sidered, can be used to specify or simplify the matter 
stress-energy tensor on £ as argued later. In summary, 
under our assumptions, the two spacetimes V ± can be 
matched within an accuracy up to the order of e 2 if any 
stress-energy tensor is accepted on £. 


E. Matching of the electromagnetic fields 


A 2 = const., 
M 2 (\ 
~R? 


l e = 


Q^_ 

2 R 2 


M\ 

a 


sin 0 cos 0 , 


_ Cn = 4 £(i + « 


Co - Co 


3 R 


R 2 R 2 J 


, M 2 ( 4 M 2 Q 2 \ , , 

k 2 (R)R — ( —1 — + -j^r ) + C 2 — £2 , 


3 R 

1 


f MR - Q 2 SM 
<0 ~^ 


ho{R) Jrn(R) { R 3 

2 (2 MR - Q 2 ){R 2 + 4 MR - 2 Q 2 ) 


+M 


6 i? 6 


+ 


m 0 


A 2 
Ao' 


h 2 (R) = 


fds(R)L 2 
M 2 (2MR — Q 2 ) 


3 R 4 


MR-Q 2 
Irn(R)R 3 


C2 


i ?£ 2 


L 2 f dS (R )' 


(2.32) 

(2.33) 

(2.34) 

(2.35) 


(2.36) 


(2.37) 


Equation (12.331) gives us the explicit form of l e , given in 
terms of the unperturbed quantities only. Note that l e 
is given as the unique equatorially symmetric solution to 
the first-order ordinary differential equation derived from 
the first junction conditions. As shown in Appendix B, 
the constant A 2 , given in Eq. (12.321) . does not appear in 
K~1. Thus, the basic properties of the thin shell matter 
are not affected by values of the constant A 2 . For the 
sake of simplicity, in this study, we chose the value of A 2 
to be 


(2.38) 


In order to determine the spacetime structure of V~, 
we have to know the specific form of the electromagnetic 
fields there. From the vector potential A+, as already 
given in Sec. II.B, we may assume the vector potential 
A~ to be given by 

A = j-^ + e 2 (B 0 (r) + B 2 {r)P 2 (cos0)) j dt ^ 
+ aA 3 (r ) sin 2 0 dcj) + 0 (e 3 ), 


where the functions A 3 , Bq, and B 2 are solutions of the 
vacuum Maxwell equations, and their explicit forms are 
derived in Appendix A. The nonzero components of the 
electromagnetic fields tangent to £ are then given by 


^ 0 $ 


F + 

r T<3> 


F, 


e$ 


+ 2 aQsin 0 cos 0 3 


-Fie = 


R 


+ 0(e 3 ), (2.40) 


2 Q (2 M 2 + 3i?£ 2 ) sin 0 cos 0 

IF 


+ 0(e 3 )(2.41) 


2a sin 0 cos 0A 3 (i?) + 0(e 3 ), (2.42) 


^r© — Fqt 

= 3e 2 AqB 2 {R) sin 0 cos 0 + 0(e 3 ), (2.43) 

If the electromagnetic fields are assumed to be regular in 
V~, the functions A 3 and B 2 , are given by 

r — L arctanh(r/T) 

A3 = C 2 


A 2 =0. 


r 


(2.44) 
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B2 = w% ( 3Lr+ ( r2 ~ 3 l 2 )f) 

3Mr 3 (2.45) 

+ ^ (2r 3 - 3 L 2 r + 3 L 3 f dS x) , 

with C '2 and a\ being constants. The e-order matching 
condition, [F©$] = 0, leads 


C 2 


Q 

R — L arctanh(i?/L) 


(2.46) 


From the e 2 -order matching condition, [Fro] = 0, we 
obtain 


a i 


2Q(2M 2 + 3 ( 2 R) 
3 Z 


2C 2 0i(3LR + (R 2 
MZ 


3 L 2 )X) 
(2.47) 


where 


X = arctanhp Z = 2R 3 - 3 L 2 R + 3 XL 3 f dS {R ). 

(2.48) 

Since the junction conditions for the electromagnetic 
fields [F Q b] = 0 are satisfied within an accuracy up to 
the e 2 order, we see that the electromagnetic fields in 
V ± can be matched through E up to the order e 2 . From 
the second equations in Eq. m , we may evaluate the 
current vectors tangent to E by using jumps of F ±a n on 
E. The nonzero components of F ±a „ are given by 


F+ T n = 


Q 


R 2 Irn 2 


e 2 Q 


GSMR 3 + 2 M 3 R - M 2 ( Q 2 + 4 R 2 ) 


6 R? (f RN f 2 
+ 6R (3 MR - Q 2 - 2R 2 ) Co 

- 2^7M 2 {R(R - 2 M) + Q 2 } 

+ 3 R {R{2R - 3 M) + Q 2 } C 2 |f 2 


0(e 3 ), 

(2.49) 


F+* = 


aQ 


R 4 fn 


T + 0(e 3 ) 


(2.50) 


:N ■ 


F 


-T 


3B' 0 (R) + 2M 2 uiA’ 3 (R) 


3A 0 fds 2 
+ {3 B' 2 {R) - 2M 2 uj 1 A' 3 (R)} P 2 


0(e 3 ), 

(2.51) 


F“* = - 


afd.s 2 A 3 (R) 

R 2 


0(e 3 )., 


(2.52) 


where Bq is the regular solution, given by 


B 0 = a 2 + 


2L0iX 

3 Mr 


(2.53) 


with a 2 being a constant. Note that a 2 is determined by 
the continuity condition of A± = on E. By us¬ 

ing these components of F ± °, we may know the surface 
current on E, j a = — [F a n ], The current vector j a satis¬ 
fies the equation of the charge conservation, i.e., j,“ = 0. 
Thus, we calculate the total charge of the thin shell E by 


Q = 


j 3 T VMld,ed<b , 


(2.54) 


where h denotes the determinant of h a &. If we assume the 
regular solutions for the electromagnetic fields, we have 


Q = Q + 0(e 3 ) , (2.55) 

because 

F~ T n = -r {3 B’ 2 (R) - 2M 2 uj 1 A' 3 (R)} P 2 + 0(e 3 ). 

3A 0 f dS i 

(2.56) 


This fact shows that our assumption that the charge of 
the black hole does not depend on the spin parameter a 
is correct within the accuracy we consider. 


F. Second junction conditions: Stress-energy 
tensor of the thin shell 

The second junction condition says that the stress- 
energy tensor of the thin shell, S a b , is given by the jump 
of the extrinsic curvature K± b across E (see Eq. (12.51) 1. 
Nonzero components of the extrinsic curvature K b up 
to the order e are given by 


k+ 1+we), 

R 3 VJRN 

(2.57) 

' / *“ + <V), 

(2.58) 

Q 2 - MR 2 , 

K T = a ~ 7^=^r + °( e )> 

VIrnR 5 

(2.59) 

r , +T 3MR-2Q 2 . 2 , 2 . 

K + Q=a — sm-0 + 0(e 2 ), 

V JRN R 3 

(2.60) 


(2.61) 

K "=-/Jk +oie% 

(2.62) 

IS—® TS — ^ \/fdS 2\ 

R e A $ R +0{€ ), 

(2.63) 

A-T = A o (^- 2 fe_) 0 + O(e »,, 

(2.64) 


(2.65) 
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The explicit forms of the e-order extrinsic curvature are 
summarized in Appendix B because their expressions are 
quite complicated. With these expressions, we may ob¬ 
tain explicit forms of the stress-energy tensor of E. 

The matter three-velocity tangent to E, u a , and the 
total energy density measured by an observer with u a , er, 
are, respectively, defined by the timelike eigenvector and 
the corresponding eigenvalue of Sg, as follows: 

S%u h = -ou a , u a u a = -1. (2.66) 

Here, u^ = e£u a means the four-velocity of the matter 
associated with the thin shell. The stress tensor of the 
thin shell, i a b, measured by an observer with u a is defined 

by 


lab = q a c q b d S cd , (2.67) 


where q a b is the projection tensor associated with u a , de¬ 
fined by q a b = S a b + u a u b . For a perfect fluid, 7 a b is pro¬ 
portional to q a b, he., lab = pqab , where p is interpreted 
as the isotropic pressure of the fluid matter. 

From the second junction condition, let us first deter¬ 
mine the explicit form of u a . The e-order relations in 
Eq. (12.661) are obviously satisfied by assuming that 


u 


a 


a 0 


(vS'"' 0 ) +o(i) - 

V fdS ~ V fltN 
4irR 


( 2 . 68 ) 

(2.69) 


where the quantity with the subscript 0 denotes the quan¬ 
tity in the no-rotation limit. With this u a , we see that 
the matter of the thin shell can be described as a perfect 
fluid within an accuracy up to the e-order. Thus, we 
have p = po + 0(e 2 ). The isotropic pressure in the limit 
of e —^ 0 is given by 


1 / M-R 1 — 2 R 2 /L 2 \ 

Po ~~8^V7 W + VTTs )' 


(2.70) 


These relations for spherically symmetric solutions may 
be obtained from the master equations given in Ref. [25[ 
if po = 0 is assumed. The solutions given in Ref. [25| can 
be then used as the unperturbed solutions for the present 
analysis if we consider the case of po = 0 (as argued in 
the later sections). Since S" 5 © = 0, the d> component of 
Eq. (12.661) leads to 


c*$ T 1 c*3> $ 3> 

+ U(fyVj — —<JU 


(2.71) 


This equation is satisfied within an accuracy up to the 
order of e if we assume u® = ttu T , where H = O(e). 
Therefore we may obtain 


(Q 2 - MR)y/Jds + (Q 2 - 2 MR)VJ^ 

\/JrnR 2 — {R 2 — 3Mi? + 2Q 2 )\fJds ' 

= afii. (2.73) 


Since 5® = 5® = 0(e 3 ), we may obtain u e = 0(e 3 ). 
Thus, we have 

u a = n T (l,0,H) + O(e 3 ). (2.74) 


The time component of the matter velocity, u T , within 
an accuracy up to the e 2 order may be obtained by the 
normalized condition of u a , 


h ab u a u b = -1. (2.75) 


The explicit form of u T is then given by 


u 


T 



+ R 4 Irn 



MR) Co 


M 2 

+R 3 SM + - \Q 2 - 2 MR 

6 1 

+4i? 2 (Q 2 - 2 MR)fl + 2 R 6 nl} 

( M 2 

+ j R(Q 2 ~ MR)( 2 + — {Q 2 - 2 MR 


(2.76) 


+2R 2 (2MR- Q 2 )Hi 



+ 0(e 3 ). 


Once u a within an accuracy up to the e 2 order is ob¬ 
tained, the e-order quantities for the thin shell may be 
calculated. The total energy density a of the shell is 
given by 


<7 = S ab u a u b , 

= ctq + e 2 Sa + 0(e 3 ), (2.77) 


where e 2 <5cr denotes the e 2 -order perturbation of the to¬ 
tal energy density a. The explicit form of 8a is given 
in Appendix C. The stress tensor of the shell, 7 ^, may 
be calculated, and its explicit expressions of its nonzero 
components are given by 


7t =Poq T T + 0(e 3 ), 7$ =Po9$ + 0(e 3 ), 

7t = Poq% + 0(e 3 ), (2.78) 






















G. A perfect-fluid thin shell 


7 + =Po 


8nK 2 


sfhsKv? + ^—^5, 


R 2 VIrn 


SM 


T 2 £o + (L 2 - 2 R 2 ) mg 

L 2 V7dS 3 

3 M 2 R + R 3 - M (Q 2 + 3 R 2 ) 


Co 


R : W Jrn 

M 2 f (R- M) (2 R 2 - 2 MR + Q 2 ) 

8R 2 fRN | 2 R 2 ^/}'rn 

(Q 2 - 2MR + i? 4 fii) {L 2 (Q 2 - 2Mi?) + 2I? 6 fii} 

l 2 r 3 V7^ 

i? (2 M 2 + 3 Q 2 ) - M (Q 2 + 4i? 2 ) 


2R 3 (MR — Q 2 ) 2 


VTr 


:N 


V/rn 

fill 


Hi 


Here and henceforth, we assume that the shell is made 
of a perfect fluid because a perfect fluid is one of the 
simplest and most useful models for describing astro- 
physical matter. For a perfect fluid, the stress tensor 
is, as mentioned before, given by 7 a b = pq a b- Since 
u e = 0 (e 3 ), as obtained in the last subsection, we obtain 
7 © = PS' 0 + 0(e 3 ) = p = po + e 2 6p , where e 2 Sp denotes 
the e 2 -order perturbation of the pressure p. Thus, the 
condition that 7 ^ is given by 7 ^ = pq a b is equivalent to 
the condition 7 ^, = pq% (see Eq. (12.7811 . in which we may 
confirm that all the components of 7^ other than 7^, are 
given in the perfect-fluid form). This condition leads 

7$ =pq$, 

= (p 0 + e 2 Sp)( 1 + zi $ u$) + 0(e 3 ), (2.81) 

= Po + (e 2 Sp + Pou®u<s>) + 0(e 3 ), 


e 2 P 2 

8nR 2 


R 2 \/JdS (h' 2 + k' 2 ) 


_ (3 R 2 - 2 L 2 ) £ 2 + {L 2 - 2R 2 ) m 2 

L 2 VJlf 

3 M 2 R + R 3 - M{Q 2 + 3 R 2 ) 

H- q -C 2 

R 3 VTrn 

M 2 \ 2 R 3 {MR-Q 2 )n 2 

~ 3 R 2 f RN \ 7f^ 

(Q 2 - 2 MR + R A Vl 1 ){L 2 (Q 2 - 2MR) + 2R 6 n 1 } 
L 2 R 3 VfdS 

R (2 M 2 + 3 Q 2 ) — M [Q 2 + 4i? 2 ) 

V7rn 

(M - R) (2 R 2 - 2 MR + Q 2 ) 

R 2 VIrn 




+ 0(e 3 ), (2.79) 


7 = 


e 2 sin 2 0 
lflTTi ? 2 


M 2 


R 2 V7rn 


M ~R 


+ Hi {M ( Q 2 + 4R 2 ) - R (2 M 2 + 3 Q 2 )} 
+ 2 R 3 ( Q 2 - MR) H 2 }- 3 ^ 2 


M 2 


VfdS 

(Q 2 - 2MR + R 4 n 1 ) 


R 3 R 2 fRN y/ fdS 

x {L 2 ( Q 2 - 2MR) + 2i? 6 Hi} 


+ 0(e 3 ), (2.80) 


where ^ = ( 7 ® ± 7q)/2. The explicit forms of ^ can 
be chosen freely as long as the first junction conditions 
for h ab given in Sec. II D are not violated. 


where the facts, u® = 0(e) and = O(e), have been 
used. Since 7 g, = Po+ c 2 Sp + 0(e 3 ), the condition (12.811) 
may be reduced to the condition for 7 ^, — 7 q , 

7%-7 % = Pou^u® + 0(e 3 ). (2.82) 

This condition and Eq. (12.801) give us £2 in terms of 
the unperturbed and first-order perturbation quantities, 
given by 


£2 = 


M 2 


^{ M -R + 2R 3 ( Q 2 - MR) fi? 

VfRN l 


3R 2 /rn 

+ (4MR 2 - 2 M 2 R + MQ 2 - 3 Q 2 R) Hi 

- (Q 2 - 2MR+ i? 4 Hi) 

L 2 ( Q 2 - 2 MR) + 2R 6 fl 1 


L 2 R 3 

+ 8nR 2 ^/JTs (Q 2 - 2.MR + i? 4 H x ) H x p 0 


(2.83) 


This condition is a third constraint for the three con¬ 
stants £ 2 , £ 2 , and D 2 besides Eqs. (12.351) and (12.371) . We 
may therefore determine values of £ 2 , £ 2 , and D 2 by solv¬ 
ing the three coupled linear algebraic equations (12.351) . 
(12.371) . and (12.831) . For a perfect fluid, the pressure is 
basically given by a function of the total energy density, 
i.e., p = p(a). The perturbations of a and p are, as seen 
from Eqs. (12.791) and (12.801) . expanded by 

Sit = 5a 0 + Sa^Pi , Sp = Sp 0 + SP 2 P 2 , (2.84) 

where Sao, Sa 2 , Spo, and Sp 2 are independent of 0. If the 
matter is described by a perfect fluid, Spo and Sp 2 have 
to relate with Sag and Scj 2 by 

(5p 0 ,5p 2 ) = (Sa 0 ,Sa 2 ) ■ (2.85) 
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As argued before, the first junction conditions and the 
assumption of the perfect-fluid shell fully determine val¬ 
ues of £ 2 , £ 2 , and D 2 . This fact directly means that St 72 
and Sp 2 are also determined. Thus, we have 


f dp\ = Sp2 

\d<7 ) 0 ( 5(72 


( 2 . 86 ) 


unless 5(i2 = 0. This equation gives us information about 
the equation of state for the thin shell matter around the 
unperturbed state given by po = p(ao). Note that in 
the present study, the two specific spacetimes are first 
matched across a thin shell E, and then the properties 
of the matter of the thin shell are determined. Thus, the 
matter properties cannot be specified before the match¬ 
ing is done. Equation (12.861) gives us a third constraint 
for SM, Co, Co, and C 4 , given by 


Spo = ~jT~ 

0(7 2 


(2.87) 


To determine the set of the constants SM, Co, Co, and 
C 4 , we need one more constraint equation for the four 
constants besides Eqs. (12.3411 . (12.361) . and (12.871) . Usually, 
we require (i) <5 <To = s, where s is a constant specified ap¬ 
propriately or (ii) the conservation of the total particle 
number of the shell as the spin parameter a increases 
from a = 0. For the requirement (i), <5 <to = s becomes 
the fourth constraint equation. By solving the four cou¬ 
pled linear algebraic equations (12.341) . (12.361) . (12.871) and 
5<Jo = s, we may obtain the set of the four constants SM, 
Co, Co, and C 4 . The simplest example is to set s = 0 . Here 
and henceforth, as unperturbed solutions, we assume the 
spherically symmetric charged regular black hole solu¬ 
tions obtained by Uchikata et al. [25|, in which the dust 
shell is considered, i.e., po = 0. The total particle number 
of the thin shell, N, is defined by 


N = 


j nu T V—hdQ d4>, 


( 2 . 88 ) 


where n is the number surface density of the fluid. The 
change in the total particle number, SN, can be given by, 

SN = J(Snu T V~-h + nSu T V—h + nu T SV~h) d<3 d<f>. 

(2.89) 

From the first law of thermodynamics, given by 


dn 

da 


a + p 


(2.90) 


and the fact that po = 0, we may obtain Sn in terms of 
5a, 


Sn = — 6a . 

<Jq 


(2.91) 


Thus, the ratio of SN to N may be explicitly given by 


SN 

~N 


Sa 0 2 Co M 2 

~a^ + ~R + TR2 
1 + Hi(2Q 2 - AMR + i? 4 Hi) 
3 R 2 f RN 


M 2 


.(2.92) 


Note that in Eq. (12.89|l . the terms related to quadrupole 
perturbations vanish after the angular integration. By 
assuming the total particle number conservation, SN = 
0, we may obtain set of solutions SM, Co, Co, and C 4 . 
The case of the requirement (ii) is investigated n detail 
in the next section. 

As argued before, the slowly rotating solutions of the 
regular black hole considered in this study are obtained 
by solving the two sets of simple coupled linear alge¬ 
braic equations whose coefficients are given by elemen¬ 
tary functions. In this study, therefore, no special tech¬ 
nique is required for the numerical calculations, all the 
numerical procedures are straightforward, and no diffi¬ 
culty appears in the numerical procedure. While any nu¬ 
merical programming languages are employable for the 
present numerical calculations, the numerical code used 
in this study is written in C. 


III. NUMERICAL RESULTS 

A. Unperturbed solutions: Spherically symmetric 
charged regular black holes 


First, we briefly describe the unperturbed regular black 
hole considered in this study. As discussed in Ref. 
[ 25 I ] . the equation of motion for a spherically symmet¬ 
ric charged dust thin shell is given by 


R+1 -§- 


q M O 2 

R + 1-—b -jp = 47rR(7o, (3.1) 


where R = dR/dr and r is the proper time of the shell. 
Note that in the present paper, M denotes the black hole 
mass, while it is used for denoting the rest mass of the 
shell in Ref. [25| . To investigate motions of the thin shell, 
it is useful to transform Eq. (13.11) to the form 

R 2 + V{R) = —1, (3.2) 


where the effective potential V is defined by 


V{R) 


+ M 
y 8nR 2 ao 


2nRaa 




R 2 

L 2 ' 


(3.3) 


Then, the stationary solution of the thin shell may be 
obtained by solving the two algebraic equations, given 

by 


V(R) = - 1, 


dV(R) 

dR 


= 0 . 


(3.4) 


simultaneously. Solutions obtained are spherically sym¬ 
metric charged regular black holes with dust thin shells 
if the radius of the shell, R, satisfies the inequalities, 
0 < R < L and R < r_, where r_ denotes the inner- 
horizon radius of the Reissner-Nordstrom black hole. For 
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4jiR 2 a 0 

FIG. 1. Radius of the unperturbed thin shell, R, given as 
a function of the rest mass of the unperturbed thin shell, 
47rf? 2 cro. Each line corresponds to the equilibrium sequence 
characterized by the same black hole mass, M, whose values 
are given near the corresponding line. 

the static solution, the derivative of the effective poten¬ 
tial, dV ( R ) / dR, can be rewritten as 

dV(R) 2 

dR ~ 

x {^^v / i^CR)+(V») ■ 

(3.5) 

Therefore, we may confirm that the condition of 
dV ( R ) /dR = 0 is equivalent to the condition of po = 0 
unless v 7 Jrn{R) = \Jfds(R)- In Figs. □ and El respec¬ 
tively, we show the radius of the shell, R, and the ratio 
of the charge to the gravitational mass of the regular 
black hole, Q/M , as functions of the rest mass of the 
shell, 47ri? 2 (7oj for the static solutions. In these figures, 
we rescale all the physical quantities by using the length 
scale of the de Sitter horizon radius L, and each line cor¬ 
responds to a sequence of the static solutions of the same 
black hole mass, M, whose values are given near the cor¬ 
responding lines. Although stable and unstable solutions 
are obtained in Ref. [2EJ, in this study, we consider the 
stable solutions only- i.e., the solutions characterized by 
d 2 V(R)/dR 2 > 0, because unstable solutions are not re¬ 
alized in nature. 

B. Dipole and quadrupole perturbations 

We numerically evaluate the perturbation quantities 
for the unperturbed regular black hole solutions, given 
in Figs. [□ and El Here and henceforth, all the perturba¬ 
tion quantities are, again, rescaled by using the length 
scale of the de Sitter horizon radius of the unperturbed 
solution, L. As shown in Sec. II.F, the dipole perturba¬ 
tions, which are of the e-order quantity and associated 
with Pi (cos d), are uniquely determined by the unper¬ 
turbed quantities only. If the stress tensor describing the 



FIG. 2. Same as Fig. [T] but for the ratio of the charge to the 
mass ratio of the black hole, Q/M. 



FIG. 3. Dimensionless angular velocity of the thin shell, 
n/e = fI/(ey/M/A 3 ), given as a function of the rest mass 
of the unperturbed thin shell, 4nR 2 ao- Each line corresponds 
to the equilibrium sequence characterized by the same black 
hole mass, M, whose values are given near the corresponding 
line. 


thin shell matter is given by the isotropic pressure only, 
as argued in Sec. II.G, the quadrupole perturbations, 
which are of the e 2 -order quantity and associated with 
P 2 (cos 0 ), are uniquely determined by the unperturbed 
and the e-order quantities. In other words, we do not 
have any degree of freedom for obtaining the dipole and 
quadrupole perturbation quantities of the regular black 
holes with perfect-fluid thin shells. To obtain the spheri¬ 
cally symmetric perturbation quantities, which are of the 
e 2 -order quantity and associated with P 0 ( cos ^)j on the 
other hand, we have options for selecting what sequences 
of slowly rotating regular black holes are considered (see 
the arguments in Sec. II.G). Therefore, we first focus on 
the results associated with the dipole and quadrupole 
perturbations. The results associated with the spheri¬ 
cally symmetric perturbations are then given in the later 
subsection. 

In Figs. 3 and 4, respectively, show the dim ensionless 
angular velocities of the shell, S2/(ey/Af/I? 3 ), and the 
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FIG. 4. Same as Fig. [3] but for the dimensionless angular 
velocity of the frame dragging in V~, w/fl. 



4itR 2 o 0 


FIG. 5. Same as Fig. [3] but for the integration constant C 2 . 


FIG. 6. Dimensionless quadrupole radial displacement of the 
thin shell defined in V~, £2 /R, given as a function of the 
rest mass of the unperturbed thin shell, 4irR 2 cro. Each line 
corresponds to the equilibrium sequence characterized by the 
same black hole mass, M, whose values are given near the 
corresponding line. 



frame dragging in V~, w/f2, along equilibrium sequences 
of unperturbed regular black holes given in Figs. [T] and 
[2j When the black hole mass and spin are kept constant, 
we observe that the dimensionless angular velocity of the 
shell, f2, increases as the radius of the shell, R, decreases 
or the charge-to-mass ratio of the black hole, Q/M , in¬ 
creases. Values of C 2 , which determine the amplitude of 
the component of the vector potential, A ^, are given in 
Fig. 5. The constant C '2 also determines the amplitude 
of the magnetic fields, B at the center of the regular 
black hole, given by 

y/WB;{r = 0) = ^\aC 2 \. (3.6) 

As for the quadrupole perturbation quantities, in Figs. 
[G] through [lT| the dimensionless radial displacements of 
the thin shell, £ 2 / R and £2 /R, the eccentricity of the thin 
shell, e/e, the dimensionless energy density perturbation 
of the thin shell, 8cr 2 /(M / A-kR 2 ), the dimensionless pres¬ 
sure perturbation of the thin shell, Sp 2 /(M/AttR 2 ), and 
the dimensionless integral constant D 2 /M 2 for the metric 
perturbations are, respectively, given as functions of the 
rest mass of the unperturbed thin shell, 47ri? 2 <70; along 


FIG. 7. Same as Fig. [6J but for the dimensionless quadrupole 
radial displacement of the thin shell defined in V + , £2 / R- 

the equilibrium sequences given in Figs. [1] and [21 Follow¬ 
ing the standard definition, see, e.g., Refs. j3jj[34j], here, 
the eccentricity of the thin shell, e, is, in terms of the 
coordinate system of V~, given by 



where r e and r p are, respectively, the effective equatorial 
and polar radii of the thin shell, given by 

r e =R + e 2 ^£0 - & , (3.8) 

r p = R -\- 1 2 (£0 + £2 + k 2 R). (3-9) 

In Figs. [6] and [7J we observe that £2 and £2 show differ¬ 
ent behaviors; values of 1/2 are always negative but the 
signs of values of £2 depend on unperturbed solutions. 
However, these behaviors are not directly related to the 
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4tiR 2 cTq 

FIG. 8. Same as Fig. ® but for the eccentricity of the thin 
shell, e/e. 



FIG. 9. Same as Fig.® but for the dimensionless quadrupole 
energy density perturbation, Sa 2 = Sat/(M/4nR 2 ). 


intrinsic physical properties of the thin shell because they 
are coordinate-dependent quantities defined in the differ¬ 
ent coordinate systems. Thus, we need quantities defined 
intrinsically in order to see intrinsic physical properties. 
To see the degree of the deformation of the thin shell, 
for instance, its eccentricity e, defined in Eq. CHI) , is a 
reasonable physical quantity. In Fig. [HI we see that the 
eccentricity of the thin shell, e, is always a real number 
and all the thin shells are oblately deformed due to ro¬ 
tation like standard rotating objects. Comparing Fig. [9] 
to Fig. [TUI we see that the equilibrium sequence charac¬ 
terized by M = 1.9 shows some peculiar behaviors. In 
these figures, we observe that near the right end point of 
the M = 1.9 sequence, 5 ct 2 becomes positive, but 6 p 2 is 
always negative along the equilibrium sequence. This is 
an unusual situation because v 2 < 0 if we naively assume 
the sound speed v s to be given by v 2 = dp/5a. Although 
this definition for the sound speed is reasonable, we can¬ 
not know how to define the actual sound speed of the 
thin shell matter at the level of the present approxima¬ 
tion because we consider stationary states only. Thus, 
the appearance of this peculiar behavior does not imme¬ 
diately mean a flaw in our treatment. Aside from the 



FIG. 10. Same as Fig. [6l but for the dimensionless quadrupole 
pressure perturbation, Sp 2 = Sp 2 / (M / An R 2 ). 



4jiR 2 o 0 


FIG. 11. Same as Fig.® but for the dimensionless integration 
constant D 2 /AF 2 . 


problem with the sound speed, the present scheme for 
obtaining the slowly rotating regular black hole with the 
perfect-fluid thin shell becomes unavailable for the unper¬ 
turbed solution having 5a 2 = 0 as discussed in Sec. II.G. 
Thus, the unperturbed solutions having 6 a 2 > 0 are not 
considered in the later discussions. For other sequences 
considered in this study, that peculiar behavior with 5 a 2 
is not observed and values of <5 02 and 5 p 2 are always neg¬ 
ative as shown in Figs. IT and ITUl The integral constant 
U 2 , given in Fig. [U] basically determines the amplitude 
of the metric perturbations in V~ (for details, see Ap¬ 
pendix A). As shown in this figure, the maximum value 
of D 2 /M 2 for unperturbed solutions considered in this 
study is given by D 2 /M 2 ~ 65 for the unperturbed solu¬ 
tion characterized by M = 1.3 and ao ~ 0, which might 
give a slightly stronger restriction for suitable values of e 
due to its slightly large maximum value of D 2 /M 2 . 
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FIG. 12. Dimensionless spherically symmetric radial displace¬ 
ment of the thin shell defined in V + , Co /R, given as functions 
of the rest mass of the unperturbed thin shell, 47ri? 2 <Jo. Each 
line corresponds to the equilibrium sequence characterized by 
the same black hole mass, M, whose values are given near the 
corresponding line. 



4jiR 2 oo 

FIG. 13. Same as Fig. 1121 but for the dimensionless spheri¬ 
cally symmetric radial displacement of the thin shell defined 
in V~, Co/A. 


C. Spherically symmetric perturbations for the 
equilibrium sequences of slowly rotating charged 
regular black holes characterized by the fixed total 
particle number of the shell 

In order to obtain spherically symmetric perturbations 
due to rotation, as mentioned before, we assume that the 
total particle number or total rest mass of the charged 
thin shell is conserved as its rotation increases from zero, 
i.e., the condition SN/N = 0. We then obtain SM, Co, 
Co and C 4 for the equilibrium sequences given in Figs. [I] 
and [2] For the M = 1.9 equilibrium sequence, as men¬ 
tioned before, (dp/da )0 cannot be defined for the unper¬ 
turbed solution with 47 rT? 2 cro ~ 0.007 because 5oi van¬ 
ishes. As a result, we cannot obtain the spherically sym¬ 
metric perturbation for the unperturbed solutions having 
47 ri? 2 < 7 o ~ 0.007. For the M = 1.9 equilibrium sequence, 
therefore, we only give results for the unperturbed solu¬ 


FIG. 14. Same as Fig. 1121 but for the dimensionless 
spherically symmetric energy density perturbation, 5ao = 
hcro/(M/47ri? 2 ). The short unlabeled line at the left top cor¬ 
responds to the result for the M=1.9 sequence. 



FIG. 15. Same as Fig. 1121 but for the dimensionless spheri¬ 
cally symmetric pressure perturbation, 5p 0 = Spo / (M/47r R 2 ). 


tions with 5o-i < —0.04 in this paper. The values of (o/R, 
Co/7?, 8<7q/(M / AnR 2 ), Spo/ (M/AttR 2 ), 5M/M, and C 4 
are given as functions of the rest mass of the unperturbed 
thin shell, 47 ri? 2 CTo, in Figs. fl2l through [171 respectively. 
As shown in Figs. fi~2l and [T3l behaviors of Co are similar 
to those of Co- This is expected by Eq. (12.34II . given by 
the first junction condition. Since the Schwarzschild-like 
coordinate system is used in V~, we may evaluate the 
mean radius of the thin shell deformed spheroidally, F, 
by r = i?(l + e 2 Co /R)- As shown in Fig. [131 Co’s always 
take positive values for the equilibrium sequences charac¬ 
terized by 6 N/N = 0. This is the usual situation because 
the centrifugal force tends to increase the equatorial ra¬ 
dius of the rotating shell and its mean radius therefore 
increases if the rest mass of the shell is kept constant. In 
Figs. UTl and fl5l we observe that values of 5a 0 and 5po are 
always negative and approach zero in the do —> 0 limit 
along the equilibrium sequences characterized by a fixed 
gravitational mass M as long as we may take the oo —■> 0 
limit. As shown in Fig. 1161 the relative changes in the 
gravitational mass of the black hole, 5M/M are always 
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FIG. 16. Same as Fig. 1121 but for the relative changes in the 
gravitational mass of the back hole, 5M/M. 



4tiR 2 (Tq 

FIG. 17. Same as Fie. 1121 but for the integration constant Ca 
appearing in ho. 


positive for the unperturbed solutions considered in this 
study. This is because the rotational energy is added in 
the thin shell but the rest mass and charge of the thin 
shell is conserved. The integration constant C 4 , given in 
Fig. [T71 basically gives a value of the metric perturbation 
ho at the center of the regular black hole. Similar to the 
values of D 2 /M 2 , slightly large values of C 4 might give 
a slightly stronger restriction for suitable values of e. 


IV. SLOWLY ROTATING REGULAR BLACK 
HOLES WHOSE THIN SHELL DISAPPEARS IN 
THE NO-ROTATION LIMIT 


The physical quantities associated with the ao = 0 cases 
correspond to those shown on the left vertical axis of all 
the figures given in Sec. III. However, we reanalyze this 
case in this section because the unperturbed solutions 
and some perturbation quantities are given in simple an¬ 
alytic forms. As discussed later, the condition for de¬ 
termining the spherically symmetric perturbation differs 
from that supposed in Sec. Ill, because the condition of 
ao = 0 is exactly assumed in this section. 

The spherically symmetric solution for this case is in¬ 
vestigated in detail by Lemos and Zanchin [24|. As ar¬ 
gued in Ref. [24j, the charge, Q, and gravitational mass, 
M , of the spherically symmetric regular black hole may 
be, in terms of the matching radius, R , and the de Sitter 
horizon radius, L , given by 


Q = 


V3R 2 

L 



(4.1) 


The regular black hole solution whose matching radius, 
R , satisfies R < r_ and R < L is allowed for the range 
of R, given by 


y/3/2 < R/L < 1, (4.2) 


as argued in Refs. [H,[25j. Thus, regular black hole solu¬ 
tions considered by Lemos and Zanchin may be specified 
if values of R/L satisfying Eq. (14.21) are given. 

Using Eqs. 63 and (14.21) for the unperturbed quan¬ 
tities of the slowly rotating regular black hole, we may 
evaluate the rotational effects on the regular black hole 
solution. The dipole perturbation quantities are then 
given by 


a 



n = 0 , 


, _ vsim ) 2 

2 R/L - X 


(4.3) 

(4.4) 

(4.5) 


where X = arctanh(R/L) as introduced before. The 
matter rotation angular velocity, f 2 , vanishes at the e- 
order approximation because there is no matter field at 
the surface of the de Sitter sphere inside the inner horizon 
for the unperturbed solution. 

For the quadrupole perturbations, we obtain 


4 R 5 f dS 
3 L 4 ’ 


(4.6) 


Finally, let us consider the case where the thin shell in¬ 
side the inner horizon disappears in the no-rotation limit, 
i.e., the case of ao = 0 and po = 0. In other words, we 
consider the case where although the de Sitter spacetime 
and the black hole spacetime are matched at some radius, 
R , there is no delta-function like matter distribution be¬ 
tween them in the unperturbed regular black hole solu¬ 
tion. Since these solutions are characterized by ao = 0, 
some parts of the results are already given in Sec. III. 


C2 = l^r R ‘lx) 2 (~ lr3 ( 9l4 - 3l2r2 + 10i?4 ) 

+ R 2 {27L 6 - 6 L 4 R 2 + L 2 R 4 + 10 R 6 )X 
- LR(27L 6 - 3 L 4 R 2 - 27L 2 R 4 + 11 R 6 )X 2 
+ 9 L 4 f 2 s (L 2 + R 2 ) 2 X 3 ) 

x (R(—3L 4 + 8 L 2 R 2 + 4 R 4 ) + 3 L 3 (L 2 - 3 R 2 )X )~ 4 , 

(4.7) 
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FIG. 18. The dimensionless integration constant, D 2 /M 2 , the FIG. 19. Same as Fig. 1181 but for the relative change in the 

eccentricity of the thin shell, e/e, the dimensionless energy gravitational mass, 5M/M, the dimensionless energy density 

density perturbation, Sa 2 = 5 a 2 /(M/ 4 nR 2 ), and the diinen- perturbation, <5<to = 8ao/(M/4nR 2 ), the dimensionless pres- 

sionless pressure perturbation, 8 p 2 = Sp 2 /(M/4-7T R 2 ), given sure perturbation, 5po = 5po/(M/4ivR 2 ), and the integral 

as functions of the unperturbed matching radius, R. constant C 4 . 




02 ° { 3 

+ R(7L 4 + 2 L 2 R 2 + 15R 4 )X - 4L 3 f dS (2L 2 - 3 R 2 )X 2 ^ 
- 8I? 4 (C 2 - 6)|(-3 L 2 R + 5 R 3 + 3 L 3 f 2 dS X)~ l . 

(4.8) 


By using the relations given in Eqs. (14.511 through (|4.8jl . 
we may obtain values of 8 a 2 and Sp 2 through the rela¬ 
tions, given by 


. 1 [ 2 R 3 f dS C 2 + 2 6 ~m 2 

2 'R'J-fdS l 3L 4 4A 2 




0.88 0.9 0.92 0.94 0.96 0.98 

R 


FIG. 20. Same as Fig. 1181 but for the dimensionless radial 
displacement of the matching radius, io/R, Co/R, £2 /R, and 
C 2 /R. 


Sp2 


R 3 / L 2 (L 2 - 2 R 2 )m 2 L 4 Sh s (h ' 2 + k' 2 ) 

2 ^ 4 v 4i?3 

(£ 4 - 6 L 2 r : 2 + 6 E 4 )C 2 + L 2 (2L 2 - 3R 2 )& 

4 R 5 Vf ^ 3 


2 L 2 - 3 R 2 \ 

3 L 2 J7Ts J ' 


(4.10) 


Although S£ = 0 in the no-rotation limit, the matter 
field in general appears on the surface of the de Sitter 
sphere as the rotation parameter, a, increases from zero. 
More specifically speaking, the energy density perturba¬ 
tion 8 a and the pressure perturbation 8 p basically appear 
at the e 2 -order approximation in the present treatment, 
as shown in Eqs. (14.91) and (14.101) . 

Values of the quadrupole perturbation quantities, the 
dimensionless integration constant, D 2 /M 2 , the eccen¬ 
tricity of the thin shell, e/e, the dimensionless energy 
density perturbation, 8a 2 /(M/4nR 2 ), and the dimen¬ 
sionless pressure perturbation, Sp 2 /(M/ 47ri? 2 ), are given 


as functions of the unperturbed matching radius, R 1 in 
Fig. M All the physical quantities are rescaled by using 
the length scale of the de Sitter horizon L. In Fig. [lHJ 
we observe that values of 8 a 2 become positive for un¬ 
perturbed solutions having 0.995 < R/L < 1 but values 
of Sp 2 are always negative. This kind of peculiar be¬ 
havior is already found in the cases of the unperturbed 
solutions with do / 0 as discussed in Sec. III. Thus, 
we henceforth focus on the unperturbed solutions with 
V3/2 < R/L < 0.995 in order to avoid problems related 
to Sa 2 = 0. 

When the spherically symmetric perturbations are con¬ 
sidered, we confront another problem with 8 a in this case 
because (To = 0. As shown in Fig. |TH1 8 a 2 s basically take 
negative values. Thus, it is possible for the total energy 
density to have a negative value somewhere regardless of 
values of e if values of 8 ao are smaller than some critical 
values. Since the appearance of negative energy density 
is unphysical, we do not consider the solutions with a < 0 
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in this study. To focus on positive energy density models 
only, in this study, we assume that 5a$ = |<5cr 2 1 if <5cr 2 < 0 
and Sag = 5(72/2 if 5a 2 > 0. Thanks to this assumption, 
Sa = 5a 0 + 5 o2P -2 (cos 6) > 0 is always satisfied. The nu¬ 
merical results of the spherically symmetric perturbation 
quantities, the relative change in the gravitational mass, 
5M/M , the dimensionless energy density perturbation, 
< 5 <To/(M/ 47 ^ii ,2 ), the dimensionless pressure perturbation, 
5pq/(M/4ttR 2 ), and the integral constant C 4 are shown 
in Fig. [19] As mentioned before, we only show the results 
for the unperturbed solutions with 0.866 < R/L < 0.99 
since values of 5o2 become positive for 0.99 < R/L. In 
this figure, we see that all the quantities given in the fig¬ 
ure take positive values except for the case of C 4 . Values 
of C 4 diverge in the limit of R/L —> 1 due to the terms 
proportional to 1 /fds, even though the results are not 
given in Fig. 1191 The dimensionless radial displacement 
of the thin shell, £,o/R, Co/A!, &/R, and C 2 /A, is given 
as functions of the matching radius of the unperturbed 
solution, R , in Fig. (2D1 


V. CONCLUSION 

We have obtained rotating solutions of regular black 
holes which are constructed of de Sitter spacetime with 
the axisymmetric stationary perturbation within the 
timelike charged thin shell and the Kerr-Newman geom¬ 
etry with sufficiently small rotation outside the shell. To 
treat the slowly rotating thin shell, we basically employ 
the method developed by de la Cruz and Israel (26| . The 
thin shell is assumed to be composed of a dust in the zero- 
rotation limit and located inside the inner horizon of the 
black hole solution. We expand the perturbation in pow¬ 
ers of the rotation parameter of the Kerr-Newman metric, 
a, up to the second order. It is found that the thin shell 
in general has anisotropic pressure at the second-order 
approximation of the rotation parameter i.e., the slowly 
rotating thin shell may not be in general composed of a 
dust with the present treatment. We may, however, set 
the shell to be composed of a perfect fluid with isotropic 
pressure by using the degree of freedom appearing in the 
physically acceptable matching of the two distinct space- 
times. 

By assuming the thin shell matter to be described by 
the perfect fluid, we investigate rotational effects on the 
regular black holes in detail. The dipole and quadrupole 
perturbations are uniquely determined by physical quan¬ 
tities of unperturbed spherically symmetric regular black 
holes and the rotation parameter, a/M. To obtain the 
spherically symmetric perturbation, we need to specify 
what equilibrium sequences of slowly rotating regular 
black holes are considered. In this study, we assume 
the following two cases: (i) Values of the spherically 
symmetric energy density perturbation are specified in 
some appropriate way. (ii) The total particle number of 
the thin shell is conserved as the rotation parameter in¬ 
creases from zero. Under the assumption (ii), we numeri¬ 


cally evaluate the perturbation quantities due to rotation 
for several equilibrium sequences of the unperturbed so¬ 
lutions of the spherically symmetric regular black hole 
characterized by a fixed gravitational mass, M, consid¬ 
ered in Ref. [25|. For the case of the equilibrium sequence 
with M = 1.9, values of < 5 <T 2 become positive for the un¬ 
perturbed solutions with 47 tI? 2 cto // 0.007, while values 
of 5p2 are always negative. This is an unusual situation, 
because it implies that dp/da < 0. The removal of this 
peculiar property with the matter field is in principle 
difficult with the present treatment, because the space- 
time metric is a prior given then the matter properties 
read out through Einstein equations. Fortunately, this 
peculiar behavior is not found for the other equilibrium 
sequences considered in this study. The present analysis 
is also applied to the case where the slowly rotating thin 
shell disappears or the stress energy tensor of the thin 
shell vanishes in the zero-rotation limit. This case cor¬ 
responds to choosing the spherically s ymm etric regular 
black hole solution considered in Ref. [24] as an unper¬ 
turbed solution. The assumption (i) is used in this case 
in order to avoid the situation where a thin shell having 
negative total energy density is obtained. 

The present results show that the regular black hole 
considered in Refs. [13, HH can rotate within an accu¬ 
racy up to the second order of the spin parameter of 
the black hole. As mentioned in Sec. Ill, the coefficients 
of the second-order perturbations in V~ appear to be 
slightly larger than those of standard rotating objects 
like rotating stars (see, e.g., Refs. [31,, HU])- This fact 
might give a slightly stronger restriction for suitable val¬ 
ues of the dimensionless spin parameter e. Therefore, it 
is interesting and important to see whether the present 
e 2 -order approximation results can be generalized to in¬ 
clude higher-order rotational effects. Studies on rotating 
regular black holes with nonperturbative approaches are 
also interesting. In this study, we focus on the case of 
the perfect-fluid thin shell. However, the thin shell with 
anisotropic pressure might be more reasonable than the 
case of the perfect fluid, because the elemental field, like 
a scalar field naturally gives the stress energy tensor with 
the anisotropic pressure. Thus, it is not unimportant to 
investigate in detail the case of the non-perfect-fluid thin 
shell. These investigations remain as future work. 
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Appendix A: Metric and Electromagnetic fields in 

V~ 

In this appendix, we briefly describe the metric 
and electromagnetic fields assumed in V~. Following 
the standard prescription given in Ref. [fiTT (see also 
Refs. (321434^ ). we consider stationary and axisymmet- 
ric perturbations around the de Sitter solution given in 
the static coordinates. The metric and electromagnetic 
field perturbations are given as solutions of the perturbed 
vacuum Einstein-Maxwell equations including the cosmo¬ 
logical constant A, 


a[c;+aj;i = a 


2 ( F^F va - i 5^F a pF a0 


A [V Q i**“] = 0, 


(Al) 

(A2) 


where G p v and V Q mean the Einstein tensor and covari¬ 
ant derivative with respect to the perturbed metric, re¬ 
spectively, and A q denotes the perturbation of the quan¬ 
tity q. In this study, as mentioned in the main text, we 
assume the equatorial symmetry and the time-azimuth 
reflection symmetry for the spacetime and electromag¬ 
netic field considered, which are reflected by the symme¬ 
try of the Kerr-Newman solution assumed as V + . Taking 
account of the rotational effects up to the second order 
of e = a/M , we may therefore assume the metric and 
vector potential as follows. 

The line element: 

ds 2 = 

- fds{r) {l + 2e 2 (ho(r) + h 2 (r)P 2 (cos 9))} dt 2 

+/„-»{ (A3) 

+ r 2 (l + 2e 2 k 2 (r)P 2 (cos 9)) 

x jd $ 2 + sin 2 9 (d<j> — 0(r)dt) 2 j + 0 (e 3 ), 

where fds(r) = 1 — r 2 /L 2 , w(r) = 0 (e), and Pi denotes 
the Legendre polynomial of degree l. Here, we employ the 
gauge condition for the metric perturbation the same as 
that of Ref. [31]. For the regular black hole considered 
in this study, we assume that 0 < r < L. 

The vector potential: 

i = |-^ +e 2 (B 0 (r) +B 2 (r)P 2 (cos 9)) J dt 
+ aA 3 (r) sin 2 9 d(j> + 0(e 3 ), 

where R and Q denote the radius and charge of the thin 
shell in the limit of e —> 0 , respectively. 

The <j> component of the Maxwell equation, 
A [V q F^“] = 0, leads the e-order equation, 


f' 

a // I J dS a / “ A _ n 

^3 + 7 ^3 7T? ^3 — u , 

JdS r fas 


(A5) 


where the prime (') means the derivative with respect to 
r. The general solution of this equation is given by 

= g 1 + C 2 ’--A»ctenl 1 ( r /£), 
r r 

where G\ and C 2 are integral constants. From the regu¬ 
larity of the electromagnetic fields at r = 0 , we need to 
impose C\ = 0. We then have 


a, = a 


r — Lx 


where 


x = arctanh(—). 

Lj 


(A7) 


(A 8 ) 


The other e-order equation is obtained from AfG 4 ^] = 0, 
given by 

(rV)' = 0 . (A9) 

From the regularity at r = 0, thus, we have 

0 = constant. (A10) 

Since the regular solutions within an accuracy up to 
the order of e have been obtained, we may move on to 
considerations of the next-order solutions. The t compo¬ 
nent of the Maxwell equation, A [W a F ta ] = 0, leads the 
e-order equations, given by 


B. 


// . 2i?2 6-B2 


r r 2 f dS 
4 lM 2 0\C 2 


2 1 L arctanh(r/T) 3 — 2 (r/L) 


2 > (All) 




r 3 f dS Zr*f 2 s 

2 B' 0 _ AM 2 0 1 C 2 


3L 2 f 


(A12) 


dS 


where 0 = a 0 i with 0 i being a constant of 0(1). The 
general solutions of these equations are given by 

a 0 L 2 f dS -3 L 2 r + 2 r 3 + 3 L 3 f d sx 

B 2 = -r-[- ai- 


2 r 3 


M 2 0iLC 2 
6 r 3 


{ ~L 2 fds log y — 4 :L 2 x + 6Lr } , (A13) 


a 3 2 LM 2 0iC 2 x 

B 0 =a 2 -1--- 

r or 


where 


V 


I ■ r 
L — r 


(A14) 


(A15) 


Here, ao, 01 , a 2 and a 3 are integral constants. Near the 
origin, r ~ 0 , these solutions become 


B 2 = -^ + *L + 0(d) 


B 0 =- 


03 


02 


20 iG2 \ 
3 M ) 


0(r 


(A16) 

(A17) 
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Therefore, if we choose 


a 0 = 0 , a 3 = 0 , 


we have the regular solutions, given by 

M 2 C 2 L0 1 r o n , 

B 2 = -^4—- {3 Lr + (r 2 - 3L 2 )x} 


(A18) 


3 r 3 

+ ( 2 r 3 - 3 L 2 r + 3 L 3 f dS x) , 

2 M 2 L 0 i x 


(A19) 


S 0 — a 2 


3r 


where we have used the relation 2x = log y. Note that 
this identity will be frequently used in the following al¬ 
gebraic manipulation. 

The other e“-order equations are obtained from the 
perturbed Einstein equations. As can be seen from 
Eq. ED, the right-hand sides of the Einstein equations 
are composed of quadratic forms in A 3 (r) and A' 3 (r). The 
independent set of the master equations for h 3 (r), mo(r), 
h 2 (r), m 2 (r), and k 2 (r) are summarized as follows: 


M 2 


m o ~ | 2 (^ 3) 2 + fdsr 2 (A ' 3 )-} , 


„ l-3(r/L ) 2 

= —IT72- m o 


(A20) 


dS 


r 2f2 

+ 4 {--W- 2 &W} 


m 2 = -rf dS <{ h 2 - -M 2 f dS (A' 3 ) 2 


(A21) 


(A22) 


,, , f 1 — 2{r/L) 2 

k 2 + h 2 = —— 2 - m 2 


fdS 


h 2 + AM 2 A 3 A’ 3 


rfdS 


3 r 2 


(A23) 


K = 


2(r/L) 2 - 1 ,, , 1 - 3(r/L ) 2 


fdS 


r 2 f? 


■ m 2 


dS 


+ 


rfdS 


rfdS 


M 


-^{r 2 (A' 3 ) 2 + 4fds\A 3 ) 2 } 


ko = 


1 


(A24) 


2 f 2 

r 2 fdS 

3 


PdS 


m 2 


3 — 4(r/L ) 2 
rfdS 

2 i- 

~Y^ k 2 
r fdS 


k’ 


+ 4l r2( ' 4 » )2 ^ 4/ « ( - 43)2 l - 


(A25) 


From Eqs. (IA20I) and (IA21I) . we obtain their general 
solutions, given by 


m 0 = C 3 + {4 L 2 rx - 2L 3 f dS x 2 


6 Lr 3 


- r 2 (2L + r logy)} , 


(A26) 


ho = C 4 


M 2 {C 2 ) 2 
6 L 3 r 4 f dS 


{—2r(L 4 — 2L 2 r 2 + 3r 4 ). 




+ Lf dS (L 4 - 3r 4 ) 

+r 2 L(L 2 + 3r 2 + Lrlogy)} , (A27) 

where C 3 and C 4 are integral constants. These functions 
can be written near the origin by 


mo = C 3 + 0(r 3 ) 
M 2 (C 2 ) 2 


hr, — 


3 L 2 


+ C A + 0(r 3 ). 


(A28) 

(A29) 


Regularity at r = 0 requires that m 3 —> 0 as r —> 0; we 
then need to set 


C 3 = 0. 

We then obtain the regular solutions: 

L 2 f dS x 2 


(A30) 


m 0 = — - 


M 2 (C 2 ) 2 (r 2 — 2L 2 


3 r 


Lr 


, _ M 2 {C 2 ) 2 f L 4 -3Lr 2 + 3r 4 

h ° =C * + uH;{ - L-r - * 

(L 4 - 3r 4 )fdsx 2 L 2 + 3r 2 
2r 2 2 


(A31) 


(A32) 


By using Eqs. (IA22D (IA25I) . we may obtain the second- 
order ordinary differential equation for h 2 , given by 

L 2 fds h 2 + 2 ^——r-—^fdsti 2 - ^2 (3T 4 - r 4 )h 2 


2{C 2 ) 2 M 2 | r 2 (3i - 4 _ ni 2 r 2 + 4t . 4) + 3L e f 3 sX 2 

(A33) 


3r 6 fds 

—2L 3 r(3L 2 - 7r 2 )f dS x} = 0 . 


The general solution of this equation is given by 


h 2 = -Di ^4 + D 2 


0 c 2 ) 2 m 


2 A/T 2 


1 


5 3 \ 3 Lf dS 

~ ^7 + » , X 


Sfds \L 2 r 


8 r 3 


' 2 ™[- 16t2A + 8i5&l2 
+r{2Lr (L 2 + 5r 2 ) - 5L 4 /j s log y}] 


(A34) 


where D\ and D 2 are integral constants. If r ~ 0, this 
solution becomes 


h 2 = — 


£>1 L 2 Di 


+ 0 (r 2 ). 


(A35) 
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To have a regular solution near the center, therefore, we 
have to impose 


Di = 0. 


(A36) 


Then, the regular solution h 2 at the origin is obtained. 
From this regular solution h 2 , the regular solutions m 2 
and fc- 2 , are given by 


m 2 = -fdsrh 2 + 


2 M 2 {C 2 ) 2 
3 r 3 

r 3 , / r 2 1 

/c 2 — r/ln “t - * 


(r - Lf dS xY 


(A37) 


2 L 2 '" 2 ' V 2 T 2 fds, 
M 2 (C 2 ) 2 f A 4 + 3L 2 r 2 - 2r 4 
6 r 2 L A f dS \ 


fdS 

-—{L 4 + 5L 2 r 2 ~2r 4 )x 
r 

^fdsL^ + 7L 2 r 2 + 2r 4 )x 2 


(A38) 


We finally obtain the set of regular solutions, given by 


h 2 — D 2 


1 


8/ dS U 2 r s 


2//^f \2 I ^ fdSX 


3 Lf dS 

8 r 3 


+ M 2 (C 2 ) 


3r 4 


Lr{L 2 + 5r 2 ) + (—5A 4 + 2Lr 2 — 5r 4 )a: t ( A onA 

3 ... > q !■ i ? (A39j 


12L 3 r 3 fds 


raa = Jl? {r(3i! 


- M ; 


- 5r 2 ) - 3L 3 x/j s } 

- 14L 2 r 2 


k 2 = 


D 2 


8r 3 L 2 
M 2 (C 2 ) 2 


j I 

, 2 /5r 4 + 14A 2 r 2 ^llA 4 

(C2) (-mv-* 

L 2 f% s x 2 7L 2 — 5r 2 \ 
3r 3 12rL 2 ) ’ 


{r(3L 2 + 4r 2 ) - 3A(A 2 + r 2 )a;} 
A 2 — 7r 2 


(A40) 


12Ar 3 
A 4 — 4A 2 r 2 — 3 r 4 
6 A 2 r 4 


1 


12 r 2 


(A41) 


( 2 )a^+? = 


R 3 \/7rn 


SM(M - R) 


M 2 


+ — {Q 4 - Q 2 R(M + 2R) - 2 MR 2 (M - 2R)j 


{Q 2 (3R 2 - 6 MR + 2Q 2 ) + MR 2 (3M - 2R)j 


+ 


P-2 


R 5 Y/rn 


M 2 (7MR — 5Q 2 ) 


C 2 


{3Q 2 (i ? 2 - 2 MR) 


Rf RN 

+2Q 4 + 3 M 2 R 2 - 2MR 3 } 


(B2) 


k+1 


1 j M 2 (2R 2 -4MR+Q 2 ) 


R 4 Vf^ { 3 R 

+ 2 5MR 2 + 2 (2Q 2 - 3 MR + R 2 ) Co 

2P 2 f M 2 (4Q 2 - 7MR + 2i? 2 ) 

R 4 VIrn I 3 R 

+ (2Q 2 — 3MR + R 2 ) g 2 |, 


(B3) 


(2) K +l _(2) 


M 2 (R 2 + RM~Q 2 ) „ 

=- _ , , — sm 2 0 . 

R 5 VJrn 


(B4) 


( 2 ) K + = 


_ 3 {M 2 (3Q 4 + Q 2 R{4R - 13 M) 


6R 7 YJrn 
+ 2R 2 (7M 2 - 6 MR + 2R 2 )) + 6A(2Q 4 
+ Q 2 R{3R - 8 M) + R 2 (9M 2 - 8 MA + 2A 2 ))Co 
+ 65MR 3 (2Q - 3Mi? + A 2 ) + 2 (.M 2 (3Q 4 
+ Q 2 R(7R - 13M) + A 2 (14M 2 - 15M A + 4i? 2 )) 
+ 3A(2Q 4 - Q 2 R{8M + 3A) 


+A 2 (9M 2 + AMR - 4 i? 2 ))C 2 ) A 2 } , 


(B5) 


Appendix B: The e 2 -order coefficients of the 
extrinsic curvature 


The e 2 -order coefficients of the extrinsic curvature 
associated with E are summarized in this appendix. 
Nonzero coefficients of are given for the sake of 

saving space, where ^AT^ is defined by 




(Bl) 


(2) f v- T _ 1 

T L 2 (\ffds) 3 


{$o - Wo - L 2 fd S h'o 
+ ($2 — m 2 — L 2 f^s h 2 ) P 2 } 


Mk-I+Wk-% 


2 (go + 77T.Q ) 

R 2 VT^ 


(B 6 ) 


+ 2 P 2 


( ^2g 2 + m 2 /-— , \ 

V R 2 ^JTs VfdS V ’ 

(B7) 



20 


^K~% -(s) K~t = 


3> 


36 


R 2 V7^ 


sin 2 0, (B 8 ) 


m K ~ 


L 2 R 2 (V7ds) 3 
- L 2 R 2 {f dS ) 2 K + 

+ (21? - 3 R 2 ) m 2 


- R 2 ) 6 + (2L 2 - 3 R 2 ) mo 

— (4L 2 — 5i? 2 )6 

L 2 R 2 (f dS ) 2 (h' 2 + 2k' 2 )}p 2 
(B9) 


Sa 2 


1 [ 26-ma (R 2 -3MR + 2Q 2 )( 2 

4t tR 2 \ y/JJs R 2 V7 rn 
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Appendix C: The explicit expressions for Sag and 8a 2 
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M 2 (Q 2 - 2MR + R 4 fh ) 2 
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2R 2 ( 6 M 2 - AMR + R 2 


24^f? 7 /feT 

+3 Q 2 R(R - AM) + 3Q 4 
+4i ? 4 (Q 2 - Mi?) 

—2i ? 6 (i?(i? - 3M) + 2Q 2 ) 


, (Cl) where Sa = Sao + 5a 2 P 2 . 
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